A New Type of Plasma Wakefield Accelerator Driven by Magnetowaves 
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We present a new concept for a plasma wakefield accelerator driven by magnetowaves (MPWA). 
This concept was originally proposed as a viable mechanism for the "cosmic accelerator" that would 
accelerate cosmic particles to ultra high energies in the astrophysical setting. Unlike the more famil- 
iar Plasma Wakefield Accelerator (PWFA) and the Laser Wakefield Accelerator (LWFA) where the 
drivers, the charged-particle beam and the laser, are independently existing entities, MPWA invokes 
the high-frequency and high-speed whistler mode as the driver, which is a medium wave that cannot 
exist outside of the plasma. Aside from the difference in drivers, the underlying mechanism that 
excites the plasma wakefield via the ponderomotive potential is common. Our computer simulations 
show that under appropriate conditions, the plasma wakefield maintains very high coherence and 
can sustain high-gradient acceleration over many plasma wavelengths. We suggest that in addition 
to its celestial application, the MPWA concept can also be of terrestrial utility. A proof-of-principle 
experiment on MPWA would benefit both terrestrial and celestial accelerator concepts. 

PACS numbers: 98.70.Sa, 52.40.Db, 52.25.Xz, 52.65.Rr 



Ultra high energy cosmic rays (UHECR) with energy 
beyond 10^"eV have been observed, yet its astrophysi- 
cal source and acceleration mechanism remain a mystery. 
This has been considered as one of the 11 most impor- 
tant questions in cosmology and particle astrophysics in 
the 21st century[l|. Plasma wakefield accelerators 
are known to possess two salient features: (1) The energy 
gain per unit distance does not depend on the particle's 
instantaneous energy or momentum. This is essential to 
avoid the gradual decrease of efficiency in reaching ultra 
high energies; (2) The acceleration is linear. Bending of 
the trajectory is not a necessary condition for this mecha- 
nism. This helps to minimize inherent energy loss which 
would be severe at ultra high energy. These qualities 
suggest that plasma wakefield might be a good candi- 
date for the cosmic accelerator. However, high-intensity, 
ultra-short photon or particle beam pulses that excite the 
laboratory plasma wakefields are not readily available in 
the astrophysical setting. In 2002 it was proposed @ that 
large amplitude plasma wakefields can instead be excited 
by the astrophysically more abundant plasma "magne- 
towaves", whose field components are magnetic in na- 
ture {\B\ > \E\). As astrophysical outfiows are most 
likely random, protons may be stochastically accelerated 
beyond ZeV energy by riding on such random wake- 
fields, which would explain the observed UHECR inverse- 
power-law energy spectrum. Recent particlc-in-cell sim- 
ulations have demonstrated that this MPWA mechanism 



is indeed vahd[6[. 

In this paper we suggest that in addition to its ap- 
plication to the celestial accelerator, this new way of 
plasma wakefield excitation may also be applicable to 
high acceleration-gradient terrestrial accelerators. In the 
case of the electron-beam-driven Plasma Wakefield Ac- 
celerator (PWFA) 01 and the laser-driven Laser Wakefield 
Accelerator (LWFA) [I, the drivers are physical entities 
that can independently exist outside of plasma. In con- 
trast the driver of the MPWA, the magnetowave, is a 
"medium wave" which cannot exist without the plasma 
medium. The production of ultra short and high inten- 
sity laser and electron beams requires complex and ex- 
pensive facilities. Magnetowaves on the other hand can 
be excited easily by transversely perturbing a background 
static magnetic field. This salient feature would hopefully 
be an attractive alternative to drive the plasma wake- 
field accelerator for laboratory applications. A proof-of- 
principle experiment on MPWA would not only benefit 
the terrestrial accelerator technology but also contribute 
to the understanding of celestial acceleration. 

Magnetized plasmas support a variety of wave modes 
propagating at arbitrary angles to the imposed magnetic 
field. For our purpose, we focus on wave propagation 
parallel to the external magnetic field to ensure the lin- 
ear acceleration. In this case, the electromagnetic waves 
become circularly polarized and the dispersion relation 
in the non-relativistic limit is 17| 
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FIG. 1: (a) Frequency and (b) phase velocity versus wavenum- 
ber for different magnetic field strengths. The vertical solid 
lines are the mean values of the pulse wavenumbers that were 
chosen for the PIC simulation for Ceise 1 and Case 2, and the 
shaded regions their ranges. 



where the upper (lower) signs denote the right-hand 
(left-hand) circularly polarized waves. ujip_ep ~ 
^yAire'^np/mi^e is the ion and electron plasma frequency 
and Wic.ec = eB/rm^eC the ion and electron cyclotron fre- 
quency, respectively. Each polarization has two real so- 
lutions with high and low frequency branches and both 
have a frequency cutoff which forms a forbidden gap for 
wave propagation. The right-hand polarized, low fre- 
quency solution is called the whistler wave which prop- 
agates at a phase velocity less than the speed of light. 
When the magnetic field is sufhciently strong such that 
ojec 2> Wep, the dispersion of the whistler mode becomes 
more linear over a wider range of wavenumbers with 
phase velocity approaching the speed of light (see Fig.l). 
The E and B components of the wave are now compara- 
ble in strength. In this regime the traveling wave pulses 
can maintain their shape over macroscopic distance, a 
condition desirable for plasma wakefield acceleration. 

Applying the dispersion relation for the whistler wave, 
with the ion motion neglected, we obtain the pondero- 
motive force acting on an individual electron as[8[ 
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where (C) is the amplitude of the whistler wave- 
packet, and C = z — Vgt the co-moving coordinate for 
the driving pulse. From here on, we drop all the in- 
dices for ion and electron. Combining this equation with 
the continuity equation and the Poisson equation, the 
longitudinal electric field in the plasma, i.e., the plasma 



wakefield, can be solved and it reads 
with 
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dC'£;^(C')cos[fcp(c-C')], (4) 



where i?,,,, is the maximum value of . An expres- 

sion similar to Eq. Q has been obtained for laser-induced 
wakefield in a magnetized plasmaQ. For a Gaussian 
driving pulse with E„[C,) = E'j^, exp(— C^/2cr^), it can 
be shown that behind the driving pulse, i. e., \C,\ ^ cr. 



X(C) = —kpfre" 



V4 



cos fcpC = X cos kpC,. 



(5) 



It is customary to express the plasma wakefield in terms 
of the Lorentz invariant "strength parameter" of the driv- 
ing pulse, flo = e_B,,,, /'^eCW, and the "cold wavebreaking" 
field, -E^ifa = m^cuipl e. Assuming the driving pulse fre- 
quency is centered around lo and its speed Vg fa Lo/k ^ c, 
the maximum wakefield. or the acceleration gradient^ at- 
tainable behind the driving pulse is then 



G = 



(1 — ujc/^y 



;XeE, 



wb 



(1 - LOc/i^y 



;XeE.u]b- (6) 



We see that in the non-relativistic limit, i.e., the linear 
regime of plasma perturbation, the wakefield in a mag- 
netized plasma is enhanced over the unmagnetized case 
by a factor 1/(1 — ojc/'^)^, which can be substantial when 
uj approaches ujc- 

We have conducted computer simulations to study 
the MPWA process driven by a Gaussian whistler pulse 
described above. Our simulation model integrates the 
relativistic Newton-Lorentz equations of motion in the 
self-consistent electric and magnetic fields determined 
by the solution to Maxwell's equations ITl| . The 
4-dimensional phase space {z,px,Py,Pz) is used for the 
charged particle dynamics and a uniform external mag- 
netic field, Bq, is imposed in the z-direction. In order to 
investigate the influence of magnetization on the wake- 
field excitation, we study two cases with the normal- 
ized physical parameters tOc/'-Op ~ 6 (Case 1) and 12 
(Case 2). Wc used a wavepackct with Gaussian width 
a — 80A/-\/2, where A is the cell size taken to be 
unity, and we assigned the wavenumber k = 27r/60A and 
27r/40A to Case 1 and Case 2, respectively. All other 
basic parameters arc common, i.e., we set mi/rrie ~ 2000 
and imposed a uniform background plasma with elec- 
tron collisionlcss skin depth, c/wp = 30A. This gives 
uj/ujp ~ 2.98 and Vg/c ~ oj/ck = 0.95 for Case 1 and 
uj/uip = 4.64 and Vg/c ~ uj/ck = 0.99 for Case 2. Other 
numerical parameters used were: total number of cells in 
the z-direction, = 2^'^ A ~ 546c/wp, average number 
of particles per cell was 10, and the time step cUpAt — 0.1 
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for Case 1 and 0.05 for Case 2. The fields were normal- 
ized by {l/30)E^b- 

To compare the two cases on the equal footing, we 
choose to fix the gradient G. This is done by setting 
the maximum amplitude for Case 1 at E^^ = 8.15 and, 
commensurate with the increase in uj, the maximum am- 
plitude for Case 2 at i?,,, = 20. These give the normal- 
ized vector potential oq — eE^ /nieCW = 0.09 for Case 1 
and 0.14 for Case 2. Thus the wakefield in our simula- 
tion is in the linear regime. The pulse was initialized at 
zq ~ 500A = IG.GGc/cjp. To avoid spurious effects, we 
gradually ramped up the driving pulse amplitude until 
t = lOOw-i for Case 1 and t = 200w-i for Case 2, dur- 
ing which the plasma feedback to the driving pulse was 
ignored. After these times, the driving pulse-plasma in- 
teraction was tracked self-consistently. As the dispersion 
relation in this regime is not perfectly linear, there was 
a gradual spread of the pulse width. Thus x and of 
the driving pulse decreased accordingly. As a result, the 
maximum wakefield amplitude, E^-, declined in time. We 
shall see, however, that such degradation becomes much 
milder as the ratio oiujc/^p increases. Fig. 2 is a snapshot 
of Ex and E^ in Case 1 at At = lQOuj~^ after the pulse 
was released. The initial wakefield amplitude agrees well 
with the theoretical maximum of E^ ~ 0.0059£'„,6. Fig. 3 
shows the snapshot of the same run at the later time 
At ~ BOOcjp"^. Here we see the obvious dispersion of the 
driving pulse. We note, however, that while the driving 
pulse continues to disperse, the wakefield remains coher- 
ent albeit with a much lower amplitude. 




FIG. 2: A snapshot of the plasma wakefield, (in black), 
induced by the whistler pulse, E^c (in gray) for Case 1 at the 
time step At = lOOcj^^. 

Next we repeat the exercise with Case 2. Fig. 4 shows 
the snapshot of Case 2 at At = 100w~^. Comparing with 
Case 1 at the same time step, we see that the wakefield 
amplitudes in Case 2 had barely degraded and were more 
constant than that in Case 1, which is indicative that the 
driving pulse had smaller dispersion for Case 2. Fig. 5 
shows the snapshot of Case 2 at At = 300w~^. Here 
the difference is more dramatic. Whereas the driving 
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FIG. 3: A late time snapshot for Case 1 at At — SOOtJp ^ . 
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FIG. 4: A snapshot of the plasma wakefield, Ez (in black), 
induced by the whistler pulse. Ex (in gray) for Case 2 at the 
time step At = 100a;j7^- 



pulse had already widely spread and the corresponding 
wakefield had largely diminished in Case 1 by this time 
step, that in Case 2 were still well maintained. 

In Fig[6l we show the PIC simulations driving pulse 
intensity contours in the to — k space for the two cases. 
There are superimposed with the theoretical curves for 
the whistler wave dispersion relations deduced from 
Eq.(l). We confirm that our driving pulses were indeed 
whistler waves. 

In our previous paperQ, we have confirmed via PIC 
simulations the concept of plasma wakefield excited by a 
magnetowave in a magnetized plasma. In this paper we 
further demonstrate that by raising the ujc/ujp ratio, the 
driving whistler pulse becomes less vulnerable to the dis- 
persion. This aspect is especially important for MPWA 
to be a viable mechanism for terrestrial accelerator since 
it is essential for an accelerated particle to continuously 
gain energy from the plasma wakefield so as to attain 
high energy. Different from the laser and the particle 
beam, the magnetowave is a medium wave which can- 
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FIG. 5: A late time snapshot for Case 2 at At = SOOcOp 
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not exist without the support from the plasma. MPWA 
should thus be of interest as a fundamental phenomena 
in plasma physics and an alternative approach to plasma 
Wakefield acceleration. 

Aside from the issues of dispersion of the driving pulse 
and sustainment of the plasma wakefield, we have not 
yet investigated other key accelerator physics issues in 
MPWA. It is curious whether MPWA would be more 
or less sensitive to plasma uniformity than PWFA and 
LWFA do. With regard to the phase slippage between 
the driving pulse and the accelerated particle, we see from 
our simulation results that by raising uJc/'-^p from 6 to 12, 
the phase velocity of the wakefield increases from 0.95 to 
0.99. This indicates that phase velocities can approach 
the speed of light if this ratio is further raised. One chal- 
lenge common to plasma wakefield accelerators is that, in 
order to concentrate the driving pulse energy, its trans- 
verse size is often as small as a few plasma skin-depths. 
This is known to induce large transverse wakefield which 
in turn creates problems with undesirably large betatron 
motions for the accelerated particle. It would be inter- 
esting to know if this new type of plasma wakefield driver 
using the magnctowaves can help to ameliorate this dif- 
ficulty. 

Finally, it would be extremely exciting if proof-of- 
principle experiments on MPWA can be pursued. With 
regard to the possible physical mechanism to excite the 
whistler magnetowave driving pulse for experimentation, 
it has been shown that a fast ion-acoustic wave can decay 
into a whistler wave plus an ion-acoustic wavep^. It is 
therefore conceivable that such a decay process, or con- 
versely the fusion of two ion-acoustic waves, can produce 
whistler wave. Inspired by this, one wonders whether 
a similar process can occur between a light wave and a 
whistler wave. If so, then perhaps a laser pulses may con- 
vert into a whistler wave pulse in a magnetized plasma 
under suitable conditions. 



FIG. 6: The intensity contours of the driving pulse as a func- 
tion of (oi, k) from PIC simulation. The light cone and the 
linear dispersion curves for the whistler wave with ujc/ujp = 1, 
6 and 12 are superimposed. 
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